Giant magneto- drag in graphene at charge neutrality 
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We report experimental data and theoretical analysis of Coulomb drag between two closely posi- 
tioned graphene monolayers in weak magnetic field. Close enough to the neutrality point, coexistence 
of electrons and holes in each layer leads to a dramatic increase of the drag resistivity. Away from 
charge neutrality, we observe non-zero Hall drag. The observed phenomena are explained by decou- 
pling of electric and quasiparticle currents which are orthogonal at charge neutrality. The sign of 
magneto-drag depends on the energy relaxation rate and geometry of the sample. 

PACS numbers: 72.80.Vp, 73.50.Jt, 73.22.Pr, 73.50.-h 



Recent measurements of frictional drag in graphene- 
based double-layer devices revealed the unexpected phe- 
nomenon of giant magneto-drag at the charge neutrality 
point. Applying external magnetic fields as weak as 0.1- 
0.3 T results in the reversal of the sign and a dramatic 
enhancement of the amplitude of the drag resistance. If 
the device is doped away from charge neutrality, the im- 
pact of such a weak field on the drag resistance is very 
modest. The observed effect weakens at low tempera- 
tures, hinting at the classical origin of the phenomenon. 

In this Letter we report experimental data on longitu- 
dinal and Hall drag resistivity in isolated graphene layers 
separated by a 1 nm thick boron-nitride (hBN) spacer. 
The observed effects are explained in terms of coexisting 
electron and hole liquids in each layer [H, Q . This theory 
is based on the hydrodynamic description of transport in 
graphene derived in Refsjj-l6| using the quantum kinetic 
equation framework 0, It provides a simplified de- 
scription of the magneto-drag effect while capturing the 
essentially classical physics of the phenomenon The 
effect can be traced back to the fact that the Lorentz 
force in the electronic band is opposite to that in the 
hole band, which is also the reason for the anomalously 
large Nernst effect 13, HI and vanishing Hall effect at 
charge neutrality. 

The classical mechanism behind the giant magneto- 
drag is illustrated in Fig. [1] The upper panel shows two 
infinite graphene layers at charge neutrality. The driving 
electric current ji in the active layer corresponds to the 
the counter-propagating flow of electrons and holes with 
zero total momentum due to exact electron-hole sym- 
metry (hence, in the absence of additional correlations 
there is no drag at the Dirac point [sl llil - fisl ) . In a weak 






FIG. 1: (Color online) Mechanism of magneto-drag at charge 
neutrality. Upper panel: in an infinite system quasi-particle 
currents in the two layers (denoted by Pi) flow in the same di- 
rection, leading to positive drag. Lower panel: in a thermally 
isolated system no net heat flow is possible; the quasi-particle 
currents in the two layers have opposite directions yielding 
negative drag. 

magnetic field, electrons and holes are defiected by the 
Lorentz force and drift in the same direction. The result- 
ing quasi-particle fiow. Pi , carries a non-zero net momen- 
tum in the direction perpendicular to the charge flow, ji. 
The momentum transfer by the interlayer Coulomb inter- 
action induces the quasi-particle current P2 in the same 
direction as Pi . The Lorentz forces acting on both types 
of carriers in the passive layer drive the charge flow in the 
direction opposite to ji . If the passive circuit is open, this 
current is compensated by a finite drag voltage, yielding 
a positive drag resistivity. 

This mechanism of magneto-drag at charge neutrality 
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is closely related to the anomalous Nernst effect in single- 
layer graphene Indeed, the quasi-particle cur- 
rent is proportional to the heat current at the Dirac point. 
A similar mechanism, where the role of Pi is played by 
a spin current, has been proposed in Ref. [l^ as a pos- 
sible explanation for a giant non-local magnetoresistance 
at charge neutrality. 

The above argument describes the steady state in the 
infinite system where all physical quantities are homo- 
geneous in real space. This is not the case in a rela- 
tively small mesoscopic sample. Whether a particular 
sample should be considered "small" or "large" is deter- 
mined by comparing the sample size to the typical length 
scale corresponding to the leading relaxation process. At 
high enough temperatures, the heat currents are most ef- 
ficiently relaxed by electron-phonon scattering, which we 
describe in this Letter by the phenomenological length 
scale £pb ^\^. 

In a finite system, the quasi-particle currents must van- 
ish at the boundaries. For W :$> £ph, the quasi-particle 
current and density is homogeneous in the bulk and the 
system remains effectively infinite. 

For W ^ ^ph, the currents Pi acquire a dependence 
on the coordinate y. In this case energy conservation 
dictates that P2{y) = —Pi{y). As the result the elec- 
tric charge in the passive layer tends to flow in the same 
direction as ji, see Fig. [TJ Thus, the drag is negative, 
which is the conventional sign for the Coulomb drag in a 
system with the same type of charge carriers. 

In order to test the above ideas, we perform new mea- 
surements of the drag effect in magnetic field which are 
illustrated in Fig. [2l The experiments are carried out on 
graphene double-layer structure with 1 nm hBN spacer 
and two electrostatic gates. The schematics of the exper- 
iment is shown in the inset of the panel D in Fig. [S] The 
same device was used in Ref. [l| for drag measurements in 
zero magnetic field. 

The map for the drag resistivity, p^xO^T: Vb), is shown 
in Fig.[2l panel A at T = 240 K. The main difference com- 
pared to zero field experiment reported earlier jH is large 
negative drag at the double neutrality point. A dramatic 
change in drag resistivity with the applied magnetic field 
is shown in more details in Fig. [51 panel B at a lower 
temperature, T — 160 K. To ensure same charge densi- 
ties ni and 712 in the top and bottom layers, we sweep 
both gates simultaneously along the line connecting the 
bottom left and top right corners of the map. The exper- 
iment shows a large negative peak in at the double 
neutrality point, ni = 712 = 0, as expected for a small 
sample (see above); in our device both layers have the 
width W ~ 2/im and sufhciently resistive contacts. 

In addition to the longitudinal drag resistivity we also 
measure the Hall drag resistivity, ^^^(Vt, Vb), shown in 
Fig. [21 panel C at T = 240 K as a function of the top 
and bottom gate voltages. Due to the low density of 
states in graphene and the small separation between lay- 



ers (in this experiment d « 1 nm), the relationship be- 
tween gate voltages and charge densities is rather non- 
trivial. To identify sign of charge carriers at each point 
in Fig. [21 panel C, we also measured resistivity maps for 
both layers. Since the resistance of graphene is peaked at 
charge neutrality, tracking the position of the resistivity 
maximum gives the lines which divide the map into the 
electron- and hole-doped parts. Such lines are shown in 
both maps, see Fig. [21 panels C and A. The observed Hall 
drag resistance is large when one of the layers is close to 
neutrality point and vanishes if two layers have the same 
charge densities with opposite signs (a white line running 
from the top left to bottom right corner). 

We now turn to the theoretical description of the drag 
effect. Consider first the Drude model for electrons and 
holes in the two layers: 



eEi + e[vie X B] 
-eEi - e[vii, X B] 

eE2 + e[v2e X B] 
-eE2 - e[v2h X B] 



Fie + eviJMi, (la) 

F,h + evih/Mi, (lb) 

F2e+eV2e/M2, (Ic) 

F2h + ev2h/M2. (Id) 



The velocities. Via with i = 1,2 and a = e,h, are the 
mean velocities of electrons and holes in the layer i, Ei,2 
and B are the electric and magnetic fields, e is the ele- 
mentary charge, and Mi, 2 are the carrier mobilities due 
to scattering on impurities. The electric j'i and quasi- 
particle Pi currents are related to Via by [4] 

3i = e{nieVie - n^hVih), Pi = UieVie + n.^Vih, (2) 

with the electron and hole concentrations riie(h) — 

J^deiyie) [e^^^/^-)/^ + l] Here iy{e) = 2\e\/7r{hvf is 
the density of states in graphene (disregarding the mag- 
netic field), and fj,i are the chemical potentials measured 
from the Dirac point. The total charge and quasi-particle 
densities are defined as rii = riie — riih and pi = riie + riih- 
In general, the frictional force acting on each type of 
carriers can be represented by the sum 

F,a^hY^ [-f1^njb{v,a-Vji,) + -ftjn.ji,{v,a+ Vjb)] , (3) 
Jb 

where the coefficients 7 appear in monolayer graphene 
due to the absence of Gallilcan invariance. The expres- 
sion ^ can be obtained by solving the quantum kinetic 
equation in the hydrodynamic approximation [l9| . The 
dimensionless coefficients 7,°^ and "ff^ are related to mi- 
croscopic collision rates f^-J]- 

For jii = Q, the first term in Eq. ^ simplifies to 



Fla 



h-f{P,-P2), 



(4) 



where 7 = h/Trp, with Tp ^ being the momentum relax- 
ation rate. The second term in Eq. ^ renormalizes the 
mobilities SSIl- The model ^ with the force dJ) also 
describes the case /ii = /12 3> T, where 7 = h/ p,iTp. 
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FIG. 2; (Color online) Panel A; Longitudinal drag resistivity in magnetic field as a function of the top and bottom gate voltages. 
Lines track positions of the neutrality points in top (open symbols) and bottom (closed symbols) layers and divide the map 
into four regions with distinct signs of charge carriers (e-e, h-h, e-h, h-e). Panel B: Magneto-drag for equal charge densities 
ni — n2 = n. Solid symbols represent the experimental data. The lines show the results of the proposed theory (see main 
text for details). Different colors correspond to diffrent values of the magnetic field. The electron-phonon relaxation length 
^ph was determined by fitting. Panel C: The map of Hall drag resistivity as function of top and bottom gate voltages. The 
white diagonal corresponds to vanishing Hall drag for ni = —n2. The lines with open and closed symbols indicate neutral 
states of the top and bottom layers, respectively. Panel D: The experimental data (blue squares, left axis) and theoretical 
plot (red solid line, right ajcis) for the Hall drag resistivity for equal charge densities n\ = n2 = n. The theoretical curve is 
calculated from the phenomenological model lU with the frictional force ((3)) as derived f9| from the quantum kinetic equation 
in the hydrodynamic approximation ^5,] using the dynamic RPA screening and the mean impurity scattering time r. Note the 
sign change at n ~ ±2 x lO^^cm"^. Inset: schematics of Hall drag measurements in double layer system. The charge density 
is controlled by voltages Vt and Vb applied to the top and bottom gate, respectively. 



For strongly-doped graphene, /i^ ^ T, the quasi- 
particle current and density are obsolete: Pi = ji/e and 
Pi — rii . Equations ([T]) may then be reduced to the stan- 
dard Drude model yielding the vanishing Hall drag re- 
sistivity p^y = —E2y/jix = and conventional drag, 
p^^ = E2x/jix = —^17/6^, which show negligible depen- 
dence on the magnetic field (20| . 

In contrast, at charge neutrality the quasi-particle and 
charge degrees of freedom are decoupled. The quasi- 
particle density for = is determined by the tem- 
perature as Pi ^ po = 7rT^/3(?iu)^, while the electric and 
quasi-particle currents ji and Pi become orthogonal, as 
shown in Fig. [TJ 



Rewriting Eqs. ([T]), (jj]) in terms of currents ([2]) and 
excluding Pi, we obtain the resistivity tensor. At charge 
neutrality, m = 0, the longitudinal drag resistivity is 
peaked and its value is given by the expression 

g _ ^7 B^M,M2 

" e2 1 + ft7Po(Mi -t- M2)/e' " 

which describes positive drag in an infinite system in 
agreement with the qualitative picture illustrated in 
Fig. [U upper panel. In the limit of weak interaction, 
< hev'^, the resuh jS]) can be obtained from 
the standard perturbative approach [l^ modified for 
graphene in a classical magnetic field. 
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Prom Eqs. Q and ([3|) we also calculate the Hall 
drag resistivity p^^. The frictional force ([3]) is derived 
[§| from the quantum kinetic equation in the hydrody- 
namic approximation ^ using the dynamic RPA screen- 
ing and the mean impurity scattering time r. The result- 
ing dependence of the Hall drag on the charge density, 
n = ni = 712, are shown in Fig. [2l panel D along with 
the corresponding experimental data without any fitting. 
The value of the scattering time t is determined from the 
single-layer resistivity measured in experiment and we 
have used the most plausible estimate for the value of the 
effective electron-electron interaction strength, a « 0.2, 
in graphene on hBN. The theory also predicts the van- 
ishing Hall drag for the case of oppositely doped layers, 
ni = —712, as well as the additional zero in the Hall drag 
at some finite value of the charge density, which agrees 
with the experimental data shown in Fig.[2l panel C, 

At the same time, the large negative peak in the longi- 
tudinal drag resistivity at the double neutrality point, 
see Fig. [21 panel B, suggests that the sample width, 
W 2 fim is relatively small as compared to at 
T = 160 K, see Fig. [U lower panel. To account for the 
finite width of the sample, we need to re-write the equa- 
tions ^ in terms of the currents ji and Pi and allow 
for the spatially varying quasi-particle density, pi (where 
we need to use the continuity equation). The resulting 
model for the first layer at charge neutrality reads (for 
the second layer one has to replace the index 1 with 2) 

- TiVpi + eniEi + [ji x B] = piFi + ePi/Mi,(6a) 
ep,E, + e[Pi xB] ^ n^F^+jjM^, (6b) 
VPi = -(pi - po)/Tph - (pi - P2)/tq. (6c) 

Here Ti = 2nn(2 cosh/ii/2r) is the mean quasi-particle 
energy. The continuity equation for the quasi-particle 
current (|6cp includes relaxation by the electron-hole re- 
combination [2] , with describing the energy loss from 
the system, which at high e noug h temperatures is domi- 
nated by phonon scattering [17|, and Tq^ describing the 
quasiparticle imbalance relaxation due to the interlayer 
Coulomb interaction. Note, that the interlayer inter- 
action does not lead to relaxation of the total quasi- 
particle current Pi + P2, which therefore vanishes if 
= 0. Near the Dirac point, energy and momen- 
tum relaxation rates coincide (in particular, tq rp). 
In doped graphene (/i 3> T) , the recombination rates are 
exponentially suppressed. 

The continuity equation for the electric current sim- 
ply reads Vji = 0, hence jt = {ji{y),0). Within lin- 
ear response, the density pi has to be substituted by 
the equilibrium value po in the products piF and piE. 
This way we obtain the linear system of differential equa- 
tions on the functions Piy{y), jix{y), and Pi{y). Since 
the charge current acquires the dependence on the trans- 
verse coordinate, we define = E2x/{jix), where 
(jix) = W-^ jixdy. 



M = 5m7Vs 
T = 240K 
a = 0.2 




FIG. 3: (Color online) The magnetic field dependence of the 
longitudinal drag resistivity at the neutrality point calculated 
from Eq. ([7]). The positive sign of the magneto-drag in weak 
fields corresponds to the limit W S> 2^ph, where iph ~ 1/im 
for the parameters of the plot. The magnetic field dependence 
of scattering rates is disregarded in the plot. 

Solving the above differential equations for rii = and 
Mi — M, we find the drag resistivity 



Pxx 



(0) 

Pxx 



F{X,Y) = 



e(x,y) 



1 + 2X + [MBf 
1 + 2X + (M5)2 *""''Q(^-^) ' 



W 
2V 



, (7a) 
(7b) 



v/(l + 2X + (MB)2)(l + r), (7c) 



where p'^xl = 1/eAfpo is the single-layer resistivity in the 
absence of magnetic field and £ph = MTTph/e is the 
phonon-scattering length. In the weak-field limit, MB <C 
1, and for W ^ £ph, the magneto-drag ([7]) is positive and 
we reproduce the result ([S]). In the absence of phonons, 
Tph oo, i.e. for the thermally insulated system, the 
drag at the Dirac point is always negative, in agreement 
with the qualitative picture of Fig. [TJ 

The result (O is illustrated in Fig. |31 where we plot 
p^j. as a function of magnetic field for different values of 
the sample width choosing realistic values for the model 
parameters T = 240 K, M = 5m^/Vs, and ^ph 1 pm. 

For doped graphene, we have solved the model Eqs. ^ 
numerically. The density dependence of the drag resis- 
tivity in magnetic field for T = 160 K is shown in Fig. [2l 
panel B along with the experimental data. The collapse 
of the theoretical curves calculated for different values of 
the magnetic field is an artifact of the phenomenological 
model [17|, which is most reliable near charge neutral- 
ity. The numerical solution also shows the suppression 
of the Hall drag as compared to the model ([T|) explaining 
the factor of 4 difference between the theory and experi- 
ment in Fig.[2)D. Based on these results, we predict that 
in wider samples, the giant magneto-drag at the Dirac 
point should become positive. We also speculate that 
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the magneto-drag at the double Dirac point may become 
positive in stronger fields due to the magnetic-field de- 
pendence of the scattering times tq, rp, and Tph- 

In conclusion, we have measured the longitudinal and 
Hall drag resistivity in double-layer graphene and pro- 
vided a theoretical description of the observed effects. 
The giant magneto-drag at the neutrality point appears 
due to the presence of two types of carriers (electrons and 
holes), which in weak magnetic fields experience a uni- 
directional drift orthogonal to the driving current. This 
effect is specific to the neutrality point, where non-zero 
drag appears despite the exact electron-hole symmetry. 
The present theory does not rely on the Dirac spectrum 
in graphene, but is equivalent to the microscopic theory 
[^, m at and far away from the charge neutrality thus 
capturing the essential physics of magneto-drag. For a 
more accurate description of the effect at intermediate 
densities, the microscopic theory should be formulated 
on the basis of the kinetic equation 
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Upon completion of this manuscript, we became aware 
of a related work by Song and Levitov [21 1 . 
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